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ABSTRACT 

The design of synthetic gene networks requires an 
extensive genetic toolbox to control the activities 
and levels of protein components to achieve desired 
cellular functions. Recently, a novel class of RNA- 
based control modules, which act through post- 
transcriptional processing of transcripts by directed 
RNase III (Rntlp) cleavage, were shown to provide 
predictable control over gene expression and unique 
properties for manipulating biological networks. 
Here, we increase the regulatory range of the Rntlp 
control elements, by modifying a critical region for 
enzyme binding to its hairpin substrates, the binding 
stability box (BSB). We used a high throughput, cell- 
based selection strategy to screen a BSB library for 
sequences that exhibit low fluorescence and thus 
high Rntlp processing efficiencies. Sixteen unique 
BSBs were identified that cover a range of protein 
expression levels, due to the ability of the se- 
quences to affect the hairpin cleavage rate and to 
form active cleavable complexes with Rntlp. We 
further demonstrated that the activity of synthetic 
Rntlp hairpins can be rationally programmed by 
combining the synthetic BSBs with a set of se- 
quences located within a different region of the hair- 
pin that directly modulate cleavage rates, providing 
a modular assembly strategy for this class of 
RNA-based control elements. 

INTRODUCTION 

The field of synthetic biology encompasses the engineering 
of new cellular functions through the design of synthetic 
gene networks. The tuning of protein levels is critical for 
proper functioning of integrated genetic networks. For 



example, the optimization of metabolic networks often re- 
quires the precise tuning and regulation of enzyme levels 
and activities to avoid undesired consequences associated 
with metabolic burden due to gene overexpression (1,2), 
the accumulation of toxic intermediates (3-5) and the re- 
direction of metabolic flux from pathways critical to cell 
growth and viability (6-8). Altered levels of protein com- 
ponents can be achieved by controlling transcription (9- 
14), post-transcriptional stability and translation (15-20) 
and protein stability (21,22). In addition, libraries of 
genetic control elements have been generated to increase 
the precision with which protein levels can be modulated 
(3,10,1 1,14,17,18). However, the majority of gene regula- 
tory tools developed to date function in bacterial hosts, 
such as Escherichia coli. Therefore, extending toolsets of 
genetic control elements to other cellular chassis is essen- 
tial to supporting the design of more complex, integrated 
genetic networks in those organisms. 

The budding yeast, Saccharomyces cerevisiae, is a 
relevant cellular chassis in industrial bioprocessing 
(23-28). The current genetic toolbox for S. cerevisiae 
gene regulation relies primarily on transcriptional control 
mechanisms such as inducible and constitutive promoter 
systems. Many inducible promoters depend on accurately 
controlling the level of the exogenously-applied inducer 
molecule, where intermediate expression levels are deter- 
mined through the partitioning of cells in the population 
between either being fully repressed or expressing the 
desired protein (29). While engineered variants have 
been constructed that offer more tunable responses to 
varying inducer concentrations (11,12,14), these systems 
can exhibit other undesirable properties, such as pleio- 
tropic effects of the inducer molecules, undesired effects 
of altering the natural regulatory networks associated with 
the native promoter system and the cost associated with 
the inducing molecule in scale-up processes. RNA-based 
control modules based on post-transcriptional mechan- 
isms may offer an advantage since their activities are 
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independent of the choice of promoter. Moreover, RNA- 
based controllers can be combined with transcriptional 
controllers to expand the design of integrated regulatory 
networks and thus provide more sophisticated control 
strategies. While post-transcriptional regulatory elements 
have been described in yeast (30-33), synthetic libraries of 
post-transcriptional control elements that exhibit tuned 
activities across a broad regulatory range had not been 
developed. 

We previously developed a novel class of RNA control 
modules that act through post-transcriptional cleavage by 
the 5. cerevisiae Rntlp enzyme (20). Rntlp recognizes 
RNA hairpins that contain a consensus AGNN tetraloop, 
which forms a predetermined fold that is recognized by the 
dsRNA binding domain (dsRBD) of Rntlp (34-36). RNA 
hairpins cleaved by Rntlp have three critical regions: the 
initial binding and position box (IBPB), comprising the 
tetraloop; the binding stability box (BSB), comprising 
the base-paired region immediately adjacent to the tetra- 
loop and the cleavage efficiency box (CEB), comprising 
the region containing and surrounding the cleavage site 
(36). Rntlp hairpins were inserted as genetic control 
elements within the 3'-untranslated region (UTR) of a 
transcript in order to direct cleavage to that region, 
thereby inactivating the transcript and lowering target 
protein levels. We designed an initial library based on ran- 
domization of the Rntlp substrate CEB and screened this 
library to identify a set of Rntlp control modules that 
tune expression levels through differential Rntlp process- 
ing rates (20). In contrast to other post-transcriptional 
control elements, which exhibit activities that vary sub- 
stantially with genetic context (37), the Rntlp control 
elements exhibited conserved gene regulatory activities 
across varied genetic contexts. The utility of the Rntlp 
control modules was demonstrated for achieving predict- 
able control over protein levels and manipulating biologic- 
al networks. 

Here, we examined the role of a different critical region 
of the Rntlp substrate, the BSB, on Rntlp processing 
efficiencies and thus gene regulatory activities. We 
generated a library of hairpins based on randomization 
of the BSB to identify sequences that modulated Rntlp 
binding affinity. Rigid structural constraints imposed by 
the BSB resulted in a low percentage (~0.6%) of sequence 
variants in the library that exhibit Rntlp binding activity, 
and thus required the development of a selection strategy 
based upon fluorescence-activated cell sorting (FACS) by 
enriching for cells exhibiting low fluorescence. In total, 16 
unique BSBs were identified that span an intermediate 
range of protein expression levels. In vitro characterization 
assays indicated that altered expression levels are due to 
the ability of BSBs to determine the hairpin cleavage rate 
and to form active cleavable complexes with Rntlp. The 
integration of the synthetic BSB sequences with different 
synthetic CEB modules demonstrated that the BSB se- 
quences function as modules that retain their relative 
activities under the context of different CEBs. Further 
characterization indicated that proportional deviation 
from the 'parent 1 BSB was inversely related to the 
strength of the coupled CEB. Our work establishes a set 
of BSB sequences and a previously developed set of CEB 



sequences as modular units that can be implemented com- 
binatorially to build synthetic Rntlp hairpins exhibiting 
finely-tuned processing properties and an extended range 
of gene regulatory activities. The development of such 
well-characterized and distinct activity modules is critical 
for synthetic biology applications, where they can be 
combined with other components to build more 
sophisticated control devices (38). 

MATERIALS AND METHODS 

Plasmid construction 

Standard molecular biology techniques were utilized to 
construct all plasmids (39). DNA synthesis was performed 
by Integrated DNA Technologies (Coralville, IA, USA) or 
the Protein and Nucleic Acid Facility (Stanford, CA, 
USA). All enzymes, including restriction enzymes and 
ligases, were obtained through New England Biolabs 
(Ipswich, MA, USA) unless otherwise noted. Pfu polymer- 
ases were obtained through Stratagene. Ligation products 
were electroporated with a GenePulser XCell (Bio-Rad, 
Hercules, CA, USA) into E. coli DH10B (Invitrogen, 
Carlsbad, CA, USA), where cells harboring cloned plas- 
mids were maintained in Luria-Bertani media containing 
50mg/ml ampicillin (EMD Chemicals). Clones were ini- 
tially verified through colony PCR and restriction mapp- 
ing. All cloned constructs were sequence verified by Elim 
Biopharmaceuticals (Hayward, CA, USA) or the Protein 
and Nucleic Acid Facility (Stanford, CA, USA). Plasmid 
maps are available in Supplementary Figure SI. 

The construction of the Rntlp characterization plasmid, 
pCS321, and the Rntlp expression plasmid, pRNTl, have 
been previously described (20). A screening plasmid 
(pCS1585) was constructed from pCS321 by replacing the 
GAL1-10 promoter with the endogenous TEF1 promoter 
(J.C. Liang, unpublished data). A second screening plas- 
mid (pCS1748) was constructed from pCS1585 by insert- 
ing an additional open reading frame (ORF) containing 
the yeast enhanced mCherry gene, ymCherry, flanked by a 
TEF1 promoter and a CYC1 terminator (J.C. Liang et al., 
manuscript in preparation). 

Insertion of engineered Rntlp substrates and appropri- 
ate controls into the 3'-UTR of yEGFP3 in pCS321 and 
pCS1585 was performed through either digestion with ap- 
propriate restriction endonucleases and ligation-mediated 
cloning or homologous recombination-mediated gap- 
repair during transformation into S. cerevisiae strain 
W303 (MATa, his3-ll,15 trpl-1 leu2-3 ura3-l ade2-l) 
through standard lithium acetate procedures (40). The 
Rntlp substrates were amplified for insertion with both 
techniques using the forward and reverse primers 
RntGap321_fwd (5'-ACCCATGGTATGGATGAATTG 
TACAAATAAAGCCTAGGTCTAGAGGCG) and 
RntGap321_rev2 (5'-TAAGAAATTCGCTTATTTAGA 
AGTGGCGCGCCCTCTCGAGGGCG), respectively. 
The Rntlp switch (RS) and the BSB variants (RS-B03, 
RS-B12) were amplified for insertion with both techniques 
using the forward and reverse primers RS_fwd (5'-ATGG 
TATGGATGAATTGTACAAATAAAGAGCCTAGG 
AAACAAACAAACTTGATGCCCTTGG) and RS_rev 
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(5'-AAATTCGCTTATTTAGAAGTGGCGCGCCCTC 
TCGAGTTTTTATTTTTCTTTTTACGATGCTGGTA 
TC), respectively. The oligonucleotide templates of RS, 
RS-B03 and RS-B12 for amplification are RS_temp 
(5'-AAACAAACTTGATGCCCTTGGCAGCCGGATG 
TCATGAGTCCATGGCATCTGGATACCAGCATCG 
T A A A A AG AAA A AT A A A) , RS-B03_temp (5'-AAACA 
AACTTGATGCCCTTGGCAGCCGGATGTTGAAAG 
TCTTCAGCATCTGGATACCAGCATCGTAAAAAG 
AAAAATAAA) and RS-B12_temp (5'-AAACAAACTT 
GATGCCCTTGGCAGCCGGATGTAGTGAGTCCA 
CTGCATCTGGATACCAGCATCGTAAAAAGAAAA 
AT AAA), respectively. In the case of digestion and liga- 
tion, the PCR products were digested with the unique re- 
striction sites Avrll and Xhol, which are located 3nt 
downstream of the yEGFP3 stop codon and upstream of 
the ADH1 terminator. Following construction and se- 
quence verification of the desired vectors, 100-500 ng of 
each plasmid was transformed into strain W303. In the 
case of gap-repair, 250-500 ng of the PCR product and 
100 ng of plasmid digested with Avrll and Xhol were 
transformed into the yeast strain. All yeast strains harbor- 
ing cloned plasmids were maintained on synthetic 
complete media with an uracil dropout solution and 2% 
dextrose at 30° C. 

Library-scale yeast transformation 

Yeast transformations with the binding library were per- 
formed on Rntlp substrates as previously described (41). 
The C13-based binding library was amplified with 
template BndLib_C13 (5'-AGCCTAGGTCTAGAGGC 
GCTATCGTGTCATGTNNNNAGTCNNNNGCATG 
GCATGATAGCGCCCTCGAGAGGG) and forward 
and reverse primers C13BindLibgap_fwd_prmr (5'-GTA 
TTACCCATGGTATGGATGAATTGTACAAATAAA 
GCCTAGGTCTAGAGGCGCTATC) and C13Bind 
Libgap_rev_prmr (5'-AATCATAAGAAATTCGCTTA 
TTTAGAAGTGGCGCGCCCTCTCGAGGGCGCTAT 
CA), respectively. The reaction was scaled-up to 800 ul to 
obtain roughly 40-50 ug of PCR product. Eight micro- 
grams of plasmid (either pCS1585 or pCS1748) was 
digested overnight with Avrll and Xhol in 400 ul total 
reaction volume. Two tubes of DNA were made with 
375 ul of PCR product (~20ug) and 150 ul of digested 
plasmid (~3 ug). A third tube acting as a negative 
control contained 450 ul of water and 75 ul of digested 
plasmid (~ 1.5 ug). Each tube was extracted with phenol- 
chloroform (1:1) and ethanol-precipitated into fresh tubes. 

Five hundred microliters of Tris-DTT [2.5 M DTT, 1 M 
Tris (pH 8.0)] was added to a 50-ml culture of yeast strain 
W303 that was grown in YPD to an OD 600 of 1.3-1.5 at 
30°C. Following 10-15 min of additional incubation, the 
cells were collected and washed in 25 ml of ice-cold Buffer 
E [10 mM Tris (pH 7.5), 2mM MgCl 2 , 270 mM sucrose] 
and washed again in 1 ml of Buffer E before being resus- 
pended to a final volume of 300 ul in Buffer E. Sixty 
microliters of this cell mixture was added to the negative 
control tube and 120 ul was added to the two tubes con- 
taining digested plasmid and the library. After allowing 
the precipitated DNA to resuspend, 50 ul of the negative 



control or the library suspension was transferred to a 
chilled 2-mm gap cuvette and electroporated (540 V, 
25 uF, infinite resistance, 2 mm gap). Each 120-ul tube of 
library suspension contained enough material for two elec- 
troporations. Following electroporation, the cells were re- 
suspended in 1 ml of pre-warmed YPD and added to a 
fresh 15-ml Falcon tube. The cuvette was washed a 
second time with a fresh 1-ml aliquot of YPD, which 
was added to the same Falcon tube. Library electropor- 
ations were collected in the same Falcon tube (8 ml total). 
The Falcon tubes were incubated with shaking for 1 h at 
30°C. After incubation, the cells were collected and resus- 
pended in 1 ml of synthetic complete media with an uracil 
dropout solution and 2% dextrose. The resuspension was 
added to 6 ml of fresh media to prepare for FACS. 

FACS and sorted library retransformation 

The transformed binding library was grown for 2-3 days 
in liquid culture. Following this growth period, the library 
and appropriate control cultures were collected and sus- 
pended in lx PBS with 1% BSA and either 7-amino- 
actinomycin D (7-AAD; Invitrogen) or 4',6-diamidino- 
2-phenylindole (DAPI; Invitrogen) was added as a viabil- 
ity stain. The cell suspension was passed through a 40-um 
cell strainer (BD Falcon) prior to analysis on a FACSAria 
or FACSAria II flow cytometry cell sorter (Becton 
Dickinson Immunocytometry Systems, San Jose, CA, 
USA). On the FACSAria, GFP was excited at 488 nm 
and measured with a bandpass filter of 530/30 nm. 
7-AAD was excited at 488 nm and measured with a 
bandpass filter of 695/40 nm. On the FACS Aria II, 
GFP was excited at 488 nm and measured with a splitter 
of 505 nm and bandpass filter of 525/50 nm. mCherry was 
excited at 532 nm and measured with a splitter of 600 nm 
and a bandpass filter of 610/20 nm. DAPI was excited at 
355 nm and measured with a bandpass filter of 450/50 nm. 
Detailed sorting procedures are presented in 
Supplementary Figures S2 and S3. The collected fractions 
were diluted to 100 ml in synthetic complete media with an 
uracil dropout solution and 2% dextrose and grown until 
an OD 60 o of ~1.5. The culture was continually 
back-diluted and grown in successively decreasing 
culture volume for 2 days at which time freezer stocks 
were made of the fractions. From the culture, 100 ul was 
collected and the library hairpins amplified by colony 
PCR with forward and reverse primers RntGap321_fwd 
and RntGap321_rev2, respectively. The PCR products 
representing the sorted binding library were recloned 
through a gap-repair method by transforming the DNA 
with the pCS321 plasmid in yeast strain W303. 

Rntlp substrate characterization assays 

S. cerevisiae cells harboring pCS321 -based plasmids were 
grown on synthetic complete media with an uracil dropout 
solution and the appropriate sugars (2% raffinose, 1% 
sucrose) overnight at 30°C. The cells were back-diluted 
the following morning into fresh media (4.5 ml total 
volume in test tubes and 450 ul in deep-well plates) to an 
optical density at 600 nm (OD 60 o) of 0.1 and grown again 
at 30°C. After 1 h, 0.5 ml (test tubes) or 50 ul (plates) of 
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20% galactose (2% final concentration) or water (non- 
induced control) was added to the cell cultures. The cells 
were grown for another 4.5 h before measuring the fluor- 
escence levels or collecting cells for RNA extraction. Cells 
harboring pCS1585-based and pCS1748-based plasmids 
followed the same procedure as pCS321 -based plasmids, 
except 2% dextrose was the only sugar in the media and 
no induction was required. 

Fluorescence quantification 

On the Quanta flow cytometer (Beckman Coulter, 
Fullerton, CA, USA), the distribution of GFP fluores- 
cence was measured with the following settings: 488-nm 
laser line, 525-nm bandpass filter and photomultiplier tube 
setting of 5.83 (pCS321 -based) and 4.50 (pCS1585-based). 
Data were collected under low flow rates until 10000 
viable cell counts were collected. For pCS321 -based 
plasmids, a non-induced cell population was used to set 
a gate to represent GFP-negative and GFP-positive popu- 
lations. For pCS1585-based plasmids, a plasmid harbor- 
ing the same backbone as pCS1585 but with no 
fluorescence gene was used to set the GFP-negative and 
GFP-positive gates. 

The LSRII flow cytometer (Becton Dickinson 
Immunocytometry Systems) was used to measure 
mCherry and GFP fluorescence from pmCh-Y-based 
plasmids. GFP was excited at 488 nm and measured with 
a splitter of 505 nm and a bandpass filter of 525/50 nm. 
mCherry was excited at 532 nm and measured with a 
splitter of 600 nm LP and a bandpass filter of 610/20 nm. 
DAPI was excited at 405 nm and measured with a 
bandpass filter of 450/50 nm. 

Quantification of ceiiuiar transcript levels 

Total RNA from S. cerevisiae was collected by a standard 
hot acid phenol extraction method (42) and followed by 
DNase I (New England Biolabs) treatment to remove 
residual plasmid DNA according to manufacturer's in- 
structions. cDNA was synthesized from 5 Lig of total 
RNA with gene-specific primers for yEGFP3 and ACT1 
(43) (rntlp_rtpcr_rev2 and ACTl_rtpcr_rev, respectively) 
and Superscript III Reverse Transcriptase (Invitrogen) 
according to manufacturer's instructions. The forward 
and reverse primers for vEGFP3 quantification are 
rntlp_rtpcr_fwd2 (5'-CGGTGAAGGTGAAGGTGATG 
CTACT) and rntlp_rtpcr_rev2 (5'-GCTCTGGTCTTGT 
AGTTACCGTCATCTTTG), respectively. The forward 
and reverse primers for ACT1 quantification are ACT1_ 
rtpcrjwd (5'-GGCATCATACCTTCTACAACGAAT) 
and ACTl_rtpcr_rev (5'-GGAATCCAAAACAATACC 
AGTAGTTCTA), respectively. Relative transcript levels 
were quantified in triplicate from three identical reactions 
from the cDNA samples by using an appropriate primer 
set and iQ SYBR Green Supermix (Bio-Rad) on an 
iCycler iQ qRT-PCR machine (Bio-Rad) according to 
the manufacturer's instructions. For each run, a 
standard curve was generated for yEGFP3 and a house- 
keeping gene, ACT1, using a dilution series for a control 
representing no insertion of an Rntlp substrate. Relative 
)EGFP3 and ACT1 levels were first individually 



determined for each sample and then the yEGFP3 values 
were normalized by their corresponding ACT1 values. 

In vitro transcription of Rntlp substrates 

All Rntlp substrates were PCR-amplified to include an 
upstream T7 promoter site using forward and reverse 
primers Rntlp-T7-PCR_fwd_prmr (5'-TTCTAATACGA 
CTCACTATAGGGACCTAGGAAACAAACAAAGTT 
GGGC) and Rntlp-T7-PCR_rev_prmr (5'-CTCGAGTT 
TTTATTTTTCTTTTTGCCGGGCG), respectively. 
From the PCR product, 1-2 (j.g was transcribed with T7 
RNA Polymerase (New England Biolabs) in the presence 
and absence of a-P 32 -GTP. The 25-ul reaction consisted of 
the following components: 1 x RNA Pol Reaction Buffer 
(New England Biolabs), 3 mM rATP, 3 mM rCTP, 3 mM 
rUTP, 0.3 mM rGFP, 1 ul RNaseOUT (Invitrogen), 
10 mM MgCl 2 , 2mM DTT, 1 ul T7 Polymerase and 
0.5 uCi oc-P 32 -GTP. Unincorporated nucleotides were re- 
moved from the reactions by running the samples through 
NucAway Spin Columns (Ambion, Austin, TX, USA) ac- 
cording to the manufacturer's instructions. 

Rntlp expression and purification 

The pRNTl plasmid was transformed into E. coli strain 
BL21 using the Z-competent E. coli Transformation Kit 
and Buffer Set (Zymo Research, Orange, CA, USA) ac- 
cording to manufacturer's instructions. Rntlp was col- 
lected as a protein extract as previously described (44). 
Briefly, an overnight culture of BL21 cells harboring 
pRNTl was back-diluted to an OD 600 of 0.5. Once the 
culture reached an OD 60 o of 1.1-1.4, it was induced with 
1 mM IPTG and grown for an additional 3 h. The cells 
were centrifuged at 2500 g for 12min at 4°C and the re- 
sulting cell pellet was frozen in a — 80°C freezer. After 
weighing the frozen cell pellet, the cells were resuspended 
in 4ml Ni 2+ buffer [25% (v/v) glycerol, 1 M NaCl, 30 mM 
Tris, pH 8.0] per gram of harvested cells. The resuspension 
was sonicated (Heat Systems-Ultrasonics, Inc.) twice with 
the following settings: 2 x 30 s, output control 5 and 50% 
duty cycle. Cellular debris was removed by centrifugation 
at 20 000 g for 30min at 4°C and the supernatant was 
filtered through a 0.2-um pore size Acrodisc 25 mm 
syringe filter (Pall Life Sciences, Ann Arbor, MI, USA). 

Rntlp was purified from the resulting supernatant with 
one 1-ml HisTrap HP column (GE Healthcare) on an 
AKTA fast protein liquid chromatography (FPLC) 
machine (GE Healthcare). Elution of the protein was per- 
formed with an imidazole concentration of 150mM in 
Ni 2 + buffer and the protein was collected in six 1-ml frac- 
tions. Protein purification was confirmed by analyzing an 
aliquot of each fraction on a SDS-PAGE gel (NuPAGE 
4-12% Bis-Tris Gel, Invitrogen) and protein function was 
confirmed by incubating an aliquot of each fraction with a 
control Rntlp substrate and analyzing the resulting 
cleavage products on an 8% denaturing polyacrylamide 
gel. Positive fractions were pooled and concentrated to 
<3-ml volume using a Centricon Centrifugal Filter 
Device (10 000 MWCO; Millipore) according to the manu- 
facturer's instructions. The concentrated protein was then 
injected into a Slide-A-Lyzer Dialysis Cassette (10000 
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MWCO; Pierce Biotechnology) and buffer-exchanged 
twice with Rntlp Storage Buffer [50% (v/v) glycerol, 
0.5M KC1, 30mM Tris, pH 8.0, 0.1M DTT, 0.1 M 
EDTA] at 4°C. The first buffer exchange took place for 
4h and the second buffer exchange occurred overnight. 
The purified Rntlp was stored in aliquots at — 20°C. 

In vitro Rntlp substrate cleavage assay 

Cleavage assays were performed on Rntlp substrates as 
previously described (44,45). Briefly, a 10-ul mixture of 
RNA and Rntlp were incubated at 30°C for 15min in 
Rntlp reaction buffer [30 mM Tris (pH 7.5), 150mM 
KC1, 5mM spermidine, 20 mM MgCl 2 , 0.1 mM DTT, 
0.1 mM EDTA (pH 7.5)]. RNA concentrations were 
varied from 0.2 to 2.0 uM and the Rntlp concentration 
was 2.3 uM. The cleavage reaction products were 
separated on an 8% denaturing polyacrylamide gel run 
at 35 W for 30min. Gels were transferred to filter paper 
and analyzed for relative substrate and product levels 
through phosphorimaging analysis on a FX Molecular 
Imager (Bio-Rad). The levels of cleaved RNA product 
were determined and fit to a Michaelis-Menten model 
using Prism 5 (GraphPad), where a relative V max was 
calculated and reported with the standard error 
determined by the fit of the model. 

In vitro Rntlp substrate mobility shift assay 

Mobility shift assays were performed as previously 
described (44,45). Briefly, a 10-ul mixture of RNA and 
Rntlp were incubated on ice for lOmin in Rntlp 
binding buffer [20% (v/v) glycerol, 30 mM Tris (pH 7.5), 
150mM KC1, 5mM spermidine, 0.1 mM DTT, 0.1 mM 
EDTA (pH 7.5)]. The RNA concentration in all samples 
was 200 nM and the Rntlp concentration ranged from 0 
to 1.7 uM. The binding reaction products were separated 
on a 6% native polyacrylamide gel run at 350 V until the 
samples entered the gel and then at 150 V for 2h. Gels 
were transferred to filter paper and analyzed for free 
RNA and RNA-Rntlp complex levels through 
phosphorimaging analysis on a FX Molecular Imager. 
The fraction of unbound RNA to total RNA was 
determined and fit to a modified Scatchard model using 
Prism 5, where a K 0 value was calculated and reported 
with the standard error determined by the fit of the model. 

RESULTS 

Design and selection of an Rntlp binding library to 
achieve tunable gene regulatory control 

Rntlp is an RNase III enzyme that cleaves hairpin struc- 
tures in S. cerevisiae. An Rntlp substrate can be divided 
into three critical regions: the IBPB, the BSB and the CEB 
(36) (Figure 1A). The BSB has a reported structural re- 
quirement in which the three nucleotides immediately 
below the tetraloop must form Watson-Crick base-pairs 
and the nucleotides in the fourth position must also 
base-pair, in either a Watson-Crick or wobble conform- 
ation, for optimal activity (46). Rntlp initially binds to the 
tetraloop and then cleaves the hairpin at two locations 



within the CEB: between the 14th and 15th nt upstream 
of the tetraloop and the 16th and 17th nt downstream of 
the tetraloop. Naturally-occurring Rntlp hairpins have 
been identified in numerous non-coding RNAs, where 
Rntlp plays a critical role in non-coding RNA processing 
and editing (47^19) and in transcripts, where Rntlp was 
shown to play a role in controlling gene expression 
(50-52). 

We previously developed a genetic system in which 
Rntlp-mediated cleavage was used to regulate gene ex- 
pression in yeast through the placement of Rntlp 
hairpins in the 3'-UTR of a target transcript (Figure IB) 
(20). We developed a set of synthetic Rntlp hairpins based 
on sequence modification within the CEB that exhibit a 
broad range of cleavage rates and thus gene regulatory 
activities. Previous in vivo work with Rntlp demonstrated 
its ability to regulate cellular RNA levels (50,53,54). The 
CEB library provided the first study linking the modifica- 
tion of a critical region of an Rntlp substrate to in vivo 
regulatory activity. Since the BSB and CEB are required 
elements for Rntlp binding and cleavage, respectively, we 
hypothesized that a similar library screening approach 
could be applied to generate synthetic BSBs exhibiting 
different Rntlp binding affinities. While prior in vitro 
work had established that mutations in the BSB could 
detrimentally affect in vitro binding and subsequent cleav- 
age (36), no studies have examined a large BSB sequence 
space for tuning gene regulatory activities in vivo. The 
synthetic CEB and BSB elements are anticipated to act 
as modular units such that they could be implemented 
combinatorially to build synthetic Rntlp hairpins exhibit- 
ing an extended range of gene regulatory activities and 
more finely-tuned levels of expression. 

We designed an Rntlp binding library based on 
randomizing the BSB (8nt) to generate Rntlp hairpins 
that exhibit different gene regulatory activities due to 
altered binding affinity between the hairpin and Rntlp 
(Figure 1C). One of the variables in the design of the 
binding library was the CEB to place within the stem, as 
we had previously described a set of synthetic CEBs with 
modified gene regulatory activity (20). As it was unknown 
how the Rntlp hairpins would respond to changes in the 
BSB, we selected a synthetic CEB contained in the hairpin 
that demonstrated the lowest level of gene expression 
(CI 3). All hairpins with active BSBs containing the CEB 
of CI 3 were expected to have comparable gene expression 
levels, improving their probability of being identified in 
the screen. C13 is also fully base-paired, such that its in- 
tegration into an Rntlp substrate stem results in a stable 
structure that is less susceptible to changes in flanking se- 
quences and has a greater probability of maintaining the 
desired hairpin structure (Figure ID). Flow cytometry 
analysis also indicates that CI 3 achieves the greatest 
population separation from a no hairpin control of all 
synthetic CEBs (Figure IE). This separation is represen- 
tative of that expected for the active and inactive binding 
populations in the binding library, thus, increasing the 
enrichment of the cell-based sort. We refer to the 
'parent' BSB that was used in the cleavage library as 
BOO, where synthetic Rntlp hairpins are identified by 
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Figure 1. Implementation of Rntlp hairpins as post-transcriptional genetic control elements and binding library design. (A) Consensus regions of an 
Rntlp hairpin. Color scheme is as follows: CEB, red; BSB, blue; initial binding and positioning box (IBPB), green. Black triangles indicate locations 
of cleavage sites. The clamp region is a synthetic sequence that acts to insulate and maintain the structure of the control element. The regulatory 
activity of our synthetic hairpins is a function of the modular elements in the CEB and the BSB. (B) Rntlp hairpins in the 3'-UTR of a gene of 
interest (goi) reduce protein levels through transcript destabilization by endonucleolytic cleavage. The hairpin contains two tunable modules: the BSB 
(blue box) and the CEB (red box). Sequence modifications within the BSB and CEB are expected to primarily tune Rntlp's ability to bind the hairpin 
(blue dial) and to cleave the hairpin (red dial), respectively. Rntlp processing of the hairpin results in unstable transcripts and thus reduced protein 
expression levels. Barrels represent the protein product of the goi; orange oval represents Rntlp. (C) Sequence and structure of an Rntlp binding 
library based on the randomization of 8nt in the BSB and containing the CI 3 CEB. (D) Sequences and structures of synthetic CEB sequences (20) 
used in the screening and characterization of the binding library. 'BOO' refers to the 'parent' BSB used in the cleavage library. (E) Histograms of GFP 
fluorescence for cell populations containing the Rntlp constructs based on the 'parent' BSB. Numbers above each peak represent the median GFP 
fluorescence for that sample normalized to a control sample with a construct containing no Rntlp hairpin. 



their CEB and BSB as Cxx-Bxx (i.e. C13-B00) or Axx- 
Bxx (i.e. A02-B00). 

Synthetic Rntlp substrates with altered binding 
affinities were identified through a cell-based fluorescence 
screen. The designed library (N8) has a diversity of 65 536 
different sequences. Due to the rigid structural require- 
ments of the BSB, we predicted that at most 384 sequences 
(0.586%) would be actively bound and subsequently 
cleaved by Rntlp. A manual plate-based screening 
strategy was previously used to efficiently identify active 
sequences from the cleavage library. However, due to the 
loose structural requirements for the CEB, the cleavage 
library had a substantially higher positive rate than that 
anticipated for the binding library, such that the 
plate-based screen was not feasible for the number of 



colonies that would have to be screened to identify a rea- 
sonable diversity of active BSBs. As a higher throughput 
method we employed FACS to efficiently identify and 
recover the low percentage of cells in the population 
with diminished fluorescence. 

The binding library was transformed into yeast through 
a library-scale gap-repair strategy and clones exhibiting 
strong gene regulatory activity were selected through 
FACS by gating for cells with low fluorescence levels 
(Figure 2A). We performed two different single FACS 
screens based on a single-color (pCS1585; yEGFP3) 
system and a two-color (pCS1748; ymCherry and 
yEGFP3) system (Liang, J.C. et al., manuscript in prepar- 
ation). Three fractions (A-C) were collected for the 
pCS1585 system around the expression level of C13-B00. 
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Figure 2. In vivo screening of an Rntlp binding library. (A) A high-throughput, in vivo, fluorescence-based screen for Rntlp hairpin activity. The 
library was cloned through gap-repair into yeast in two different plasmid systems. Clones exhibiting low GFP fluorescence were sorted from 
the population through FACS. A sorted library pool was generated through colony PCR from collected cellular fractions and gap-repaired into 
the characterization plasmid. Clones that maintained low GFP fluorescence levels were selected for sequencing and further characterization. 
(B) FACS procedure for the single-color (pCS1585-based) system. Fractions A, B and C were collected based on exhibiting GFP levels similar to 
the median fluorescence of C13-B00 (left panel). Only fraction A maintained a low level of GFP expression after the fractions were regrown (right 
panel). (C) FACS procedure for the two-color (pCS1748-based) system. Fraction D was collected based on a gate set to collect all cells exhibiting 
GFP fluorescence levels below a cells containing a positive GFP control construct lacking an Rntlp hairpin module (left panel). Fraction D 
maintained a low level of GFP expression when regrown (right panel). 



Following the initial sort, each fraction was regrown and 
subsequently characterized by flow cytometry (Figure 2B 
and Supplementary Figure S2). Only fraction A retained a 
low level of fluorescence, whereas the regrown populations 
in fractions B and C shifted substantially towards inactive 
(high fluorescence) levels. The results suggest that sorted 
fractions B and C contain a large percentage of false posi- 
tives, or clones harboring plasmids containing inactive 
Rntlp hairpins that exhibited low fluorescence levels due 
to noise in gene expression profiles. Based on the fraction 
profiles, fraction A was selected for further testing. Most 
of the false positives due to genetic noise were removed 
with the two-color pCS1748 system. A gate representing 
diminished GFP fluorescence was determined based on a 
GFP positive control lacking an Rntlp hairpin. All library 
clones with decreased GFP levels were collected into a 
single fraction (D). Fraction D was regrown, characterized 
by flow cytometry and demonstrated to retain low fluor- 
escence levels (Figure 2C and Supplementary Figure S3). 
Therefore, fraction D was also selected for further 
characterization. 



After the completion of the FACS screens, the sorted 
constructs were recloned to remove false positives due to 
mutations in the plasmid or the yeast background that 
would cause reduced GFP levels independent of Rntlp 
activity. We retrieved the selected Rntlp hairpin se- 
quences from fractions A and D by colony PCR and 
gap-repaired the recovered hairpin constructs into 
pCS321. Individual clones were initially characterized for 
gene regulatory activity by measuring cellular fluorescence 
through a plate reader assay. Colonies positive for GFP 
knockdown were sequenced to determine the BSB 
sequence. In total 16 unique BSB sequences were identified 
including the 'parent' BSB (Supplementary Figure S4 and 
Supplementary Table SI). The predicted secondary struc- 
ture of the hairpins was determined by RNAstructure 
(http://rna.urmc.rochester.edu/RNAstructure.html). The 
binding library structures deviate from structural require- 
ments that were previously established through in vitro 
studies (46). All of the BSB structures contain Watson- 
Crick base pairing in the first 3nt below the tetraloop; 
however, in the fourth position certain sequences exhibit 
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mismatching. The results from our cell-based BSB library 
screen suggest that the in vivo structural requirements for 
the BSB are not as stringent as previously described. 

A synthetic Rntlp binding library exhibits a range of gene 
regulatory activities in vivo 

The initial screen for active BSB sequences was performed 
in the context of a CEB exhibiting high cleavage activity 
(CI 3). As a result, most of the recovered binding library 
hairpins exhibited low gene expression levels, such that 
differences in activity between the synthetic BSBs were 
difficult to resolve with the flow cytometry assay at low 
fluorescence levels (Supplementary Table SI). To gain 
better resolution on the differences in BSB activity and 
examine the activity of the synthetic BSBs in the context 
of a different CEB, we integrated the selected BSB se- 
quences within the context of an Rntlp hairpin containing 
a different synthetic CEB (A02) (Figure ID). The range 
of regulatory activities spanned by the binding library in 
the context of the A02 CEB was measured at the protein 
and transcript levels. Flow cytometry analysis of the 
synthetic Rntlp hairpins indicated that the selected set 
of hairpins spanned an intermediate gene regulatory 
range— from 25% (A02-B05) to 75% (A02-B01) (Table 
1 and Figure 3A). The regulatory activities of the selected 
hairpins are not evenly distributed across this range, with 
the majority exhibiting activities in the range of 25^5% 
relative protein levels. The results suggest that the binding 
library achieves a smaller regulatory range than that 
observed with the cleavage library and may be more ap- 
propriate for fine-tuning gene expression (20). We built 
negative controls for several binding library hairpins 
and for the 'parent' hairpins (A02-B00 and C13-B00) 
by mutating the tetraloop sequence (CAUC or GAAA) 
to impede Rntlp activity while maintaining the secondary 
structure of the hairpins. The negative controls 



Table 1. In vivo characterization data for the binding library 



Substrate 


Normalized 


Normalized 




protein 


transcript 




levels (%) 


levels (%) 


A02-B00 


28 ± 1 


43 ± 8 


A02-B00 (CAUC) 


81 ± 2 


106 ± 11 


A02-B01 


75 ± 3 


82 ± 8 


A02-B02 


62 ± 2 


64 ± 6 


A02-B03 


50 ± 2 


53 ± 5 


A02-B04 


32 ± 1 


52 ± 1 


A02-B05 


25 ± 0 


39 ± 3 


A02-B06 


27 ± 2 


57 ± 2 


A02-B07 


37 ± 3 


51 ± 3 


A02-B08 


30 ± 2 


53 ± 4 


A02-B09 


36 ± 3 


56 ± 5 


A02-B10 


42 ± 3 


55 ± 4 


A02-B11 


32 ± 2 


51 ± 4 


A02-B12 


27 ± 2 


47 ± 5 


A02-B13 


39 ± 4 


53 ± 5 


A02-B14 


48 ± 4 


70 ± 2 


A02-B15 


48 ± 4 


58 ± 7 


No insert 


100 ± 3 


100 ± 8 



All normalized protein and transcript levels were determined as 
described in Figure 3A and B, respectively. 



demonstrated that the majority of knockdown obser- 
ved from each hairpin is due to Rntlp processing 
(Figure 3B). 

The reduced protein expression levels observed from the 
Rntlp binding library is expected to be due to a reduction 
in the steady-state transcript levels due to rapid degrad- 
ation of the transcript following endonucleolytic cleavage 
by Rntlp. We measured relative transcript levels for each 
Rntlp hairpin by quantitative real-time PCR (qRT-PCR) 
(Table 1). A plot of normalized yEGFP3 expression levels 
versus normalized yEGFP3 transcript levels indicates that 
there is a strong positive correlation (r = 0.847) between 
the two measures of activity (Figure 3C). A preservation 
of rank order was also observed between protein and tran- 
script levels as indicated by the Spearman rank correlation 
coefficient (p = 0.668). Specifically, with decreasing tran- 
script levels a similar decrease in protein levels is generally 
observed, confirming that the fluorescence observed was 
due to changes in the steady-state transcript levels. The 
negative controls based on mutating the tetraloop con- 
firmed that Rntlp cleavage is the cause of the observed 
transcript knockdown (Figure 3B). 

Synthetic BSBs exhibit modular activity with different 
CEBs in vivo 

We next examined the gene expression data for the Rntlp 
hairpins harboring the synthetic BSBs in the context of 
two CEBs (CI 3, A02) for trends in regulatory activity 
across the binding library. We defined a new variable z c 
as the ratio of the knockdown from a binding library 
member (Bxx) to that of the 'parent' BSB (BOO) for a 
specific CEB (Cxx or Axx): 

knockdown c _Bxx 
z c (Bxx) = — 

knockdown f „Boo 

A z-value greater than unity indicates increased 
knockdown due to the synthetic BSB, whereas a z-value 
less than unity indicates decreased knockdown due to the 
BSB. We calculated z-values for each BSB in the context 
of the CEBs CI 3 and A02 (z cl3 and z A02 , respectively) and 
plotted the two variables against each other (Figure 3D). 
For ease in interpretation, we divided the graph into four 
regions with the point (1,1) at the intersection of the quad- 
rants. Regions I and IV indicate BSBs for which activities 
relative to parent are conserved between the different 
CEBs, whereas Regions II and III indicate BSBs that 
exhibit varying activities in the context of different 
CEBs. Nearly all BSBs are located in Regions I and IV, 
with the majority falling in Region I. The data indicate 
that if a BSB causes increased knockdown in the context 
of one CEB, it will likely exhibit the same activity in the 
context of another CEB. 

To further examine the gene regulatory activities of the 
synthetic BSBs, we determined a new variable Az c , which 
is the difference between the z c value of the 'parent' 
(which by definition is 1) and the z-value of the BSB: 

Az c (Bxx) = z c (B00) - z c (Bxx) = 1 - z c (Bxx) 

We calculated Az for each BSB in the context of the CEBs 
CI 3 and A02 and plotted the variables against each other 
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Figure 3. In vivo characterization of the selected Rntlp binding library and demonstration of the modularity of the BSB sequences. (A) The gene 
regulatory activities of the binding library spans an intermediate range of protein expression levels. Normalized protein expression levels were 
determined by measuring the median GFP levels from a cell population containing the appropriate construct through flow cytometry analysis 
and values are reported relative to that from an identical construct lacking a hairpin module. The 'parent' BSB is indicated in gray. (B) The 
transcript and protein levels associated with several binding library members and their corresponding mutated tetraloop (CAUC or GAAA) controls 
support that the observed gene regulatory activity is due to Rntlp processing. Transcript levels were determined by measuring transcript levels of 
yEGFPS and a house-keeping gene, ACT1, through qRT-PCR and normalizing the yEGFPS levels with their corresponding ACT1 levels. 
Normalized transcript levels for each construct are reported relative to that from an identical construct lacking a hairpin module ('no insert'). 
(C) Correlation analysis of protein and transcript levels from the binding library members demonstrates a strong correlation between the two 
measures of gene regulatory activity. (D) The gene regulatory activity of synthetic BSBs is conserved in the context of different CEB modules. The 
ratio of the knockdown exhibited from a binding library hairpin to that exhibited from the 'parent' (BOO BSB), r c , was determined in the context of 
two different CEBs and plotted against each other. Regions I and IV represent hairpins whose activities relative to 'parent' remain consistent when 
combined with different CEBs. Regions II and III represent hairpins whose activity varies relative to 'parent' when combined with different CEBs. 
(E) Synthetic BSBs modules generally exhibit higher relative activities in the context of weaker CEB modules. A variable representing the departure 
from 'parent' activity, Azc, was calculated in the context of two different CEBs and plotted against each other. The solid line indicates where the 
values of Azc are the same for both hairpins. 



(Figure 3E). Data points that fall on the Az ci3 = Az A02 
line would indicate BSBs that have the same proportional 
effect on knockdown for both CEBs. The data fall beneath 
the Az cl3 = Az A0 2 line in the region where Az Au2 > Az ci3 , 
indicating that Rntlp hairpins with weaker CEBs are 
affected more by changes in binding affinity through 
modification to the BSBs. The data exhibit a strong 
positive correlation (r = 0.946) and can be fit with a 



trendline by linear regression that passes close to the 
origin (0,0), suggesting that we see a consistent ratio 
between Az cl3 and Az A0 2 values, where this ratio is de- 
pendent on the CEBs. The data also exhibit a strong pres- 
ervation of rank order (p = 0.929), demonstrating that the 
relative activity between BSB modules is maintained re- 
gardless of the CEB module present in the hairpin. Taken 
together, the results show the maintenance of BSB activity 
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in connection with different CEB stems and also that the 
proportional deviation from the 'parent' BSB is 
determined by the strength of the CEB. 

We further examined the ability of the BSB modules to 
tune activities of more complex RNA control devices 
based on regulating Rntlp processing. Previous work 
demonstrated that incorporation of an RNA sensor com- 
ponent, or aptamer, within the CEB region of an Rntlp 
hairpin can result in an Rntlp switch (38). Specifically, 
a theophylline aptamer was directly integrated into the 
CEB region of an Rntlp hairpin and ligand binding was 
demonstrated to inhibit Rntlp cleavage activity, resulting 
in ligand-responsive control over expression levels. Due 
to the placement of the aptamer in the CEB, the BSB 
modules become the only library of parts that can be 
used to tune the ability of Rntlp to bind and cleave this 
sensing-actuation device. We demonstrated the activity of 
the BSB modules within the different stem/CEB context of 
this Rntlp switch. We incorporated two BSB modules 
(B03, B12) that exhibited different activities into the 
Rntlp switch hairpin (Supplementary Figure S5A). The 
switches incorporating the BSB modules (RS-B03, 
RS-B12) exhibited altered expression levels relative to 
parent that correlated with the activities observed from 
these modules in the context of the A02 stem 
(Supplementary Figure S5B). These results further support 
that the synthetic BSBs are effective in tuning expression 
levels in a variety of Rntlp hairpin stems and devices. 

In vitro characterization demonstrates that Rntlp binding 
library members achieve differential activity through 
alterations in Rntlp cleavage rates and affinity 

We hypothesized that the variation in transcript process- 
ing and subsequent protein expression levels exhibited by 
the binding library is due to variations in binding affinity 
resulting from alterations in the BSB sequence and/or 
structure. To examine whether the synthetic Rntlp 
binding library members exhibit differences in binding af- 
finity to Rntlp, we performed in vitro binding assays with 
purified Rntlp. Binding reactions were ran with 20 nM of 
in vitro synthesized radiolabeled RNA encoding an Rntlp 
hairpin and varying concentrations of purified Rntlp. The 
reactions were ran in the absence of magnesium and other 
divalent metal ions that are essential for cleavage to allow 
Rntlp to bind to the substrates without subsequent 
cleavage (55). Bound products were separated by non- 
denaturing polyacrylamide gel electrophoresis and quanti- 
fied through phosphorimaging analysis (Figure 4A). We 
analyzed the reaction through a modified Scatchard 
equation in which the fraction of unbound RNA (R) to 
total RNA (R 0 ) is plotted against the enzyme (E) concen- 
tration. The equation is as follows: 

Z = — - Kd 
R 0 K D +[E] 

The dissociation constant, K D , for each synthetic Rntlp 
hairpin was determined through this analysis method 
(Table 2). The data indicate that there is a moderate 
positive correlation (r = 0.486) between K D and in vivo 
gene regulatory activity (Figure 4B). While we observe 



several data points demonstrating similar transcript 
levels for different K G values, hairpins that bind less 
tightly to Rntlp (i.e. higher ^ D ) generally tend to have 
higher transcript levels as anticipated. The binding library 
has an expanded range of Kd values compared to the 
cleavage library, due to several library members having 
K D values greater than those previously reported with 
the cleavage library (20). In vivo, we observe that most 
binding library members have increased gene expression 
levels greater than 'parent'. In vitro, this same phenom- 
enon is experienced as most binding library members have 
decreased affinity for Rntlp. Binding library members 
that have decreased gene expression levels than 'parent' 
in vivo also have K D values comparable or less than 
'parent' in vitro. However, we also observed that the 
mutant tetraloop control bound Rntlp with a similar 
Ky> as the library hairpins, although cleavage was not 
evident (Figure 3B). It has been previously reported that 
Rntlp is able to bind its substrates in inactive and active 
conformations in vitro (45). Therefore, it is plausible that 
the binding observed in these in vitro assays is due to both 
types of complexes with Rntlp. Under this situation, the 
reported K D may not be solely related to complexes that 
can be processed. It has also been shown that changes in 
the BSB affect both binding affinity and hairpin process- 
ing by Rntlp in vitro (36). Therefore, it is important 
to examine the effects of the binding library on Rntlp 
processing rates. 

We analyzed the cleavage reaction between Rntlp and 
the binding library through a Michaelis-Menten model, 
with the substrate (S) being the hairpin transcript, the 
enzyme (E) being Rntlp and the product (P) being the 
cleaved pieces of the transcript. Under these conditions, 
the following reaction occurs: 

E+S % ES^E+P 

The rate of product formation (V) is modeled as: 

* [S] k 2 * [E] 0 * [S] 
~ K M +[S] ~ K M +[S] 

The maximum rate of product formation (F max ) is the 
product of the total enzyme concentration ([E] 0 ) and k 2 . 
Alterations in the cleavage efficiency will have an effect on 
the value of k 2 and thus F max . We performed in vitro RNA 
cleavage reactions with a constant concentration of 
purified Rntlp against a range of in vitro synthesized 
radiolabeled Rntlp hairpins to determine the relative 
values of k 2 for each synthetic Rntlp hairpin. Reaction 
products were separated by denaturing polyacrylamide gel 
electrophoresis and quantified through phosphorimaging 
analysis (Figure 4C). The resulting data were fit to the 
Michaelis-Menten model to calculate a relative cleavage 
rate (RCR), which is directly proportional to K max . The 
RCR value for A02-B00 is set to 1 and the rest of the 
reported values are normalized to A02-B00. 

The RCR values for each synthetic Rntlp hairpin 
were determined through this analysis method (Table 2). 
The results confirm that the mutant tetraloop control is 
not processed by Rntlp in vitro (supporting in vivo 
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Figure 4. In vitro characterization of the binding library demonstrates that the observed tuning of gene regulatory activity is achieved through 
modulation of cleavage rates and binding affinities. (A) Representative mobility shift assays and analyses by non-denaturing polyacrylamide gel 
electrophoresis of two binding library members: A02-B05 and A02-B14. The top band corresponds to RNA-Rntlp complexes; the bottom band 
corresponds to unbound RNA. Rntlp was added to the following final concentrations in each reaction (left to right; inuM): 0, 0.42, 0.83, 1.25 and 
1.66. (B) Correlation analysis of binding affinity (K D ) and yEGFP3 transcript levels indicates a moderate positive correlation between binding affinity 
and gene regulatory activity. The data point for the A02-B0O (GAAA) negative control is indicated in gray. (C) Representative cleavage reaction 
assays and analyses by denaturing polyacrylamide gel electrophoresis on hairpins A02-B00 (GAAA), A02-B06 and A02-B07. The top band cor- 
responds to unreacted full-length RNA; the bottom band corresponds to the three cleavage products expected from Rntlp processing. The three 
cleavage products differ in size by 1 nt and cannot be resolved into individual bands under the conditions used for this assay. RNA was added to the 
following final concentrations in each reaction (left to right; inuM): 0.2, 0.5, 1.0 and 2.0. Reactions lacking Rntlp were performed with 0.2 uM of 
RNA. (D) Correlation analysis of relative cleavage rate (RCR) and yEGFPS transcript levels demonstrates a moderate positive correlation between 
cleavage rate and gene regulatory activity. 



observations). There is a moderate positive correlation 
(r = 0.480) between the measured RCR and gene regula- 
tory activity for the synthetic Rntlp hairpins (Figure 4D). 
Generally, increases in Rntlp's ability to cleave a sub- 
strate result in greater transcript knockdown. Compared 
to the cleavage library (20), we observe a smaller range of 
RCR values with the binding library, due to the observed 
decreased range in transcript knockdown. In fact, we 
observe that cleavage library members exhibiting gene 
regulatory activities within the range exhibited by the 
binding library members have similar RCR values. 
However, a correlation analysis between the K D and 
RCR values indicates that there is no correlation 
between the binding affinity and relative cleavage rates 
for the binding library members (data not shown). 
Changing the BSB does result in changes in binding 
affinity, but it can also result in variation in the cleavage 
rate. For any given hairpin, the two properties contribute 
to the observed transcript levels. In particular, the BSB 
may affect the hairpin's ability to form active or inactive 
complexes or the BSB may affect the hairpin's processing 
rate. 



DISCUSSION 

We utilized a cell-based library screening approach to de- 
velop a set of synthetic BSB sequences to modulate the 
gene regulatory activity of engineered Rntlp hairpins. 
Previous in vitro studies showed that the BSB contains 
nucleotides critical to Rntlp binding as mutations in the 
region resulted in reduced affinity (36). However, this 
earlier work examined a relatively small sequence space 
for the BSB region and did not link the in vitro processing 
activities to in vivo gene regulatory activities. The in vitro 
studies established a consensus BSB structure, where three 
Watson-Crick base-pairs were required immediately 
below the tetraloop followed by another base-pair that 
could include the wobble guanine-uracil pair. Based on 
these reported structural requirements, we estimated that 
a small percentage (~0.6%) of the randomized BSB 
library would contain hairpins cleavable by Rntlp. 
Therefore, a high-throughput FACS-based screen was em- 
ployed to enrich the library for hairpins resulting in 
reduced GFP fluorescence. The BSB library was designed 
in the context of a synthetic CEB from the cleavage library 
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Table 2. In vitro characterization data for the binding library 
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that produced the greatest amount of knockdown (CI 3) 
(20). This design biased the library such that any positive 
hits would exhibit the lowest possible expression levels to 
enhance the separation, and thus the selection, of the 
population of cells containing active hairpins versus the 
larger population containing inactive hairpins (~99.4% of 
the population). In total, 16 unique BSBs were identified. 
However, in contrast to previous in vitro work (36), 
several of the synthetic BSBs did not contain a base-pair 
in the fourth position from the tetraloop, suggesting that 
this structural requirement is relaxed in vivo. 

The selected BSB sequences were further characterized 
in the context of a CEB that exhibited weaker gene regu- 
latory activity (A02) (20) to better resolve differences in 
the BSB activities. The regulatory activities of the BSB 
sequences are distributed across an intermediate range of 
25-75% relative protein levels with the majority exhibiting 
activities in the range of 25^45% . Interestingly, the 
binding library exhibited a decreased range of activity 
relative to the cleavage library (20). The relative effects 
of each region on Rntlp-based gene regulatory activity 
would not have been predicted a priori from prior 
in vitro characterization studies. Our results suggest that 
the binding library may be more appropriate for the 
fine-tuning of gene expression. For example, due to its 
wide range of gene expression, the cleavage library can 
be employed initially to first identify regulatory ranges 
of interest. More focused regulatory activities can then 
be explored through implementation of the synthetic 
BSBs with appropriate synthetic CEBs, which will allow 
a more careful searching of activity space around a 
selected CEB. As such, the combinatorial application of 
the synthetic BSB and CEB elements can be used to extend 
and tune the regulatory range accessible through the en- 
gineered Rntlp hairpins. To demonstrate the ability to 
combine BSB and CEB sequences, we examined the syn- 
thetic BSB sequences in the context of two synthetic CEBs 
generated in the cleavage library. Our experimental results 
show that the synthetic CEB and BSB elements act as 



modular units, that the BSBs maintain their activity 
under the context of different CEBs, and that the propor- 
tional deviation from the 'parent' BSB is determined by 
the strength of the CEB (Figure 3D and E). 

In vitro characterization studies determined the relation- 
ship between binding affinity, cleavage rate and 
gene-regulatory activity for the binding library. Changes 
in the BSB sequence are expected to result in changes in 
Rntlp cleavage rate as well as to affect Rntlp binding 
(36). We observed a moderate correlation between the 
binding affinity and transcript levels and between the 
cleavage rate and transcript levels for the binding library 
(Figure 4B and D). In contrast, although slight changes in 
affinity were detected with the cleavage library, these 
changes were not correlated with gene regulatory activity 
(20). However, there was a stronger correlation observed 
between cleavage rate and transcript levels for the cleavage 
library. The data indicate that nucleotide modifications in 
the BSB cause changes in affinity and cleavage rate; 
however, the exact contribution of each of these variables 
to the observed gene regulatory activity is unclear. In 
addition, Rntlp is known to bind in active and inactive 
complexes where the inactive complex is magnesium- 
independent (45,56). The inactive complex may be more 
prevalent in the binding assay due to the absence of mag- 
nesium in the reaction buffer, which is also critical to the 
proper folding of RNA molecules (57). The changes in K D 
are reflective of how much total RNA is bound regardless 
of conformation. Thus, it is possible that the changes to 
the BSB are affecting the partitioning between inactive 
and active states. 

This work extends the regulatory capacity of the first set 
of post-transcriptional stability control elements in yeast 
by developing a set of BSB modules that can be integrated 
with a previously described set of CEB modules to ration- 
ally design synthetic Rntlp substrates. With 16 CEB and 
16 BSB modules developed, 256 different Rntlp hairpins 
can be generated to achieve finely-tuned levels of gene 
expression in S. cerevisiae. The ability to tune gene expres- 
sion levels has been demonstrated to be critical to the op- 
timization of encoded behaviors in complex genetic 
networks in synthetic biology (58-61) and in metabolic 
engineering, where in the latter applications such tuning 
strategies play a key role in minimizing metabolic burden 
(1-2) and accumulation of toxic intermediates (3,5) and 
redirecting metabolic flux from native pathways required 
for cell growth and viability (8). In addition, having 
well-characterized sets of modules that correspond to dif- 
ferent regions of the control element is critical to extend- 
ing this modular design strategy to more complex 
sensing-actuation devices (38). For example, the BSB 
modules become the critical tuning modules in this 
device architecture, as the sensor component is located 
within the CEB region and cannot be replaced with the 
synthetic CEB modules (Supplementary Figure S5). Thus, 
in combination with other RNA components as well as 
engineered or native promoter systems, the engineered 
Rntlp control modules provide a powerful tool for 
programming genetic regulatory networks in yeast, and 
thus advancing the application of this cellular chassis in 
biomanufacturing and biosynthesis processes. 
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